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Abstract Since the first discovery of microlensing events nearly two decades ago, 
gravitational microlensing has accumulated tens of TBytes of data and developed into 
a powerful astrophysical technique with diverse applications. The review starts with a 
theoretical overview of the field and then proceeds to discuss the scientific highlights. 
(1) Microlensing observations toward the Magellanic Clouds rule out the Milky Way 
halo being dominated by MAssive Compact Halo Objects (MACHOs). This confirms 
most dark matter is non-baryonic, consistent with other observations. (2) Microlensing 
has discovered about 20 extrasolar planets (16 published), including the first two 
Jupiter-Saturn like systems and the only "cold Neptunes" yet detected. They probe 
a different part of the parameter space and will likely provide the most stringent test 
of core accretion theory of planet formation. (3) Microlensing provides a unique way 
to measure the mass of isolated stars, including brown dwarfs to normal stars. Half 
a dozen or so stellar mass black hole candidates have also been proposed. (4) High- 
resolution, target-of-opportunity spectra of highly-magnified dwarf stars provide in- 
triguing "age" determinations which may either hint at enhanced helium enrichment 
or unusual bulge formation theories. (5) Microlensing also measured limb-darkening 
profiles for close to ten giant stars, which challenges stellar atmosphere models. (6) 
Data from surveys also provide strong constraints on the geometry and kinematics of 
the Milky Way bar (through proper motions); the latter indicates predictions from cur- 
rent models appear to be too anisotropic compared with observations. The future of 
microlensing is bright given the new capabilities of current surveys and forthcoming 
new telescope networks from the ground and from space. Some open issues in the 
field are identified and briefly discussed. 



1 INTRODUCTION 



Gravitational microlensin g in the local gro up refers to the temporal brightening of a background star 
due to intervening objects. lEinsteinld 19361) first studied (micro)lensing by a single star, and concluded 
that "there is no great chance of observing this phenomenon." Although there were some works in 
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intervening years by lRefsdal d 19641) and lLiebesI (1 19641) . the field was revitalized by lPaczyfiskil dl986l) 
who proposed it as a method to detect MAssive Compact Halo Objects (MACHOs) in the Galactic 
halo. 

From observations of microw ave background radiation and nucleosynthesis (see, e.g. 
IKomatsu et alj|20llt IStei gman 2003), it is clear that most of the dark matter must be non-baryonic, 
and so the original goal of microlensing is now obsolete. Nevertheless, microlensing has devel- 
oped into a powerful technique with diverse applications in astrophysics, including constraints on 
MACHOs, the study of the structure of the Milky Way, stellar atmospheres and the detection of 
extrasolar plan ets and stellar-mass black hole candidat es. Since the first discoveries of microlensing 
events in 1993 dAlcock et al]l993"tlUdalski et al.ll 19931) . the field has made en ormous progress in the 
last two decades. A number of review s have been written on this topic (e.g. iPaczvhski 1996 : Mao 
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Evanj2003t|Wambsgansd2006h . with the most recent highlights given in Maold2008al) . lGouldl 



d2009t) and iGaudil d20ldT The readers will also greatly benefit from two recent, comprehensive 
conference proceedings: the Manchester Microlensing Conferencf] and the 2011 Sagan Exoplanet 
Summer Workshop: Exploring Exoplanets with Microlensing. The workshop materials contain not 
only recent scientific highlights but also hands-on exercises for data reduction and modelling. 

The structure of t his review is a s follows. Sec.[2]introduces the basics of gravitational microlens- 
ing, which reproduces lMaol d2008 b) in a slightly modified form; Sec.[3]builds on the introduction and 
discusses the applications of gravitational microlensing. We finish this review with an outlook for 
the field in Sec. [4] Due to the rapid expansion of the field, it is unavoidable that the reference list is 
incomplete (and somewhat biased). 



2 BASICS OF GRAVITATIONAL MICROLENSING 
2.1 What is gravitational microlensing? 

According to General Relativity, the light from a background source is deflected, distorted and 
(de)magnified by intervening objects along the line of sight. If the lens, source and observer are suf- 
ficiently well aligned, then strong gravitational lensing can occur. Depending on the lensing object, 
strong gravitational lensing can be divided into three areas: microlensing by stars, multiple-images 
by galaxies, and giant arcs and large-separation lenses by clusters of galaxies. For microlensing, the 
lensing object is a stellar-mass compact object (e.g. normal stars, brown dwarfs or stellar remnants 
[white dwarfs, neutron stars and black holes]); the image splitting in this case is usually too small 
(of the order of a milli-arcsecond in the local group) to be resolved by ground-based telescopes, thus 
we can only observe the change in magnification as a function of time. 

The left panel in Fig. Q] illustrates the geometry of microlensing. A stellar-mass lens moves 
across the line of sight toward a background star. As the lens moves closer to the line of sight, its 
gravitational focusing increases, and the background star becomes brighter. As the source moves 
away, the star falls back to its baseline brightness. If the motions of the lens, the observer and the 
source can be approximately taken as linear, then the light curve is symmetric. Since the lensing 
probability for microlensing in the local group is of the order of 10~ 6 (see Sec. 12. 81 ), the microlens- 
ing variability usually should not repeat. Since photons of different wavelengths follow the same 
propagation path (geodesies), the light curve (for a point source) should not depend on the color. 
The characteristic symmetric shape, non-repeatability, and achromaticity can be used as criteria to 
separate microlensing from other types of variable stars (exceptions to these rules will be discussed 
in Sec.|231). 

To derive the characteristic light curve shape shown in the right panel of Fig. Q] we must look 
closely at the lens equation, and the resulting image positions and magnifications for a point source. 

1 http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=54 

2 http://nexsci.caltech.edu/workshop/201 1/ 



2 



Mao 




-2-10 1 



(t-tj/t, 



Fig. 1 The left panel shows a side-on view of the geometry of microlensing where a lens 
moves across the line of sight toward a background source. The right panel shows two light 
curves corresponding to two dimensionless impact parameters, uq = 0.1 and 0.3. The time 
on the horizontal axis is centered on the peak time to and is normalized to the Einstein 
radius crossing time t.E- The lower the value of uq, the higher the peak magnification. For 
the definitions of uq and see Sec. 12. 5. ll 



2.2 Lens equation 

The lens equation is straightforward to derive. From Figure [2] of the lensing configuration, simple 
geometry yields 

V + D ds a = £^, (1) 

where D s and D^s are the distance to the lens (deflector), distance to the source and distance 
between the lens (deflector) and the source, t] is the source position (distance perpendicular to the 
line connecting the observer and the lens), £ is the image position, and a is the deflection angle. 
For gravitational microlensing in the local group, Euclidean geometry applies and D ( \ s = D s — 
Mathematically, the lens equation provides a mapping from the source plane to the lens plane. The 
mapping is not necessarily one-to-one. 

Dividing both sides of Equation ([]]) by D s , we obtain the lens equation in angles 

/3 + a = 6, (2) 

where (3 = r)/D s , 9 = £/Da, and a = a x Da s /D s is the scaled deflection angle. These angles 
are illustrated in Fig. [2] 

For an axis-symmetric mass distribution, due to symmetry, the source, observer and image po- 
sitions must lie in the same plane, and so we can drop the vector sign, and obtain a scalar lens 
equation: 

p + a = 0. (3) 

2.3 Image positions for a point lens 

For a point lens at the origin, the deflection angle is given by 

4GM 1 ^ 
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and the value of the scaled deflection angle is 



D d D ds AGM 6 E 2 

a = -Fr\ a \ = -R — ^FT77 = ^r> £ = D d 6, (5) 



D s 1 1 D s c?D, 



where we have defined the angular Einstein radius as 



= * „ 0.55 J 1 -^' f ^ f^U . (6) 
D d Y D d /Ds \8kpcJ \0.3MqJ 

The lens equation for a point lens in angles is therefore 

(3 + °-f- = 9. (7) 

We can further simplify by normalizing all the angles by 6s, r s = (3/6-e, and r = 6 /6s, then the 
above equation becomes 

1 

r s + - = r, (8) 
r 

where r s is the source position and not to be confused with the size of the star, which we denote as 
r*. 

For the special case when the lens, source and observer are perfectly aligned (r s = 0), due to 
axis-symmetry along the line of sight, the images form a ring ("Einstein" ring) with its angular size 
given by Equation ([6j. 

For any other source position r s ^ 0, there are always two images, whose positions are given by 



r s ± ^r s 2 + 4 

r± = . (9) 

The '+' image is outside the Einstein radius (r + > 1) on the same side as the source, while the '— ' 
image is on the opposite side and inside the Einstein radius (r_ < and |r_| < 1). The angular 
separation between the two images is 



A6 = 6 E (r+ - r_) = feV4 + r s 2 . (10) 
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The image separation is of the same order as the angular Einstein diameter when r s <> 1, and thus 
will be in general too small to be observable given the typical seeing from the ground (~ one arc- 
second); we can only observe lensing effects through magnification. One exception may be the VLT 
interferometer (VLTI) which can potentially resolve the two images. This may be important for dis- 
covering stellar-mass black holes since t hey have larger image separations due to their large r masses 
than typical lenses with mass ~ 0.3M Q dDelplancke eTZIlioOlHRattenburv & Maokood) . 
The physical size of the Einstein radius in the lens plane is given by 

r E = D d 9 E = « 2 . 2 Au . / 4 x °* ( 1 - °±) ( ^ ( -JL.) ^ . 

(11) 

So the size of the Einstein ring is roughly the scale of the solar system, which is a coincidence that 
helps the discovery of extrasolar planets around lenses. 



2.4 Image magnifications 

Since gravitational lensing conserves surface brightnesfl the magnification of an image is simply 
given by the ratio of the image area and source area, which in general is given by the determinant of 
the Jacobian in the lens mapping (see Sec. 12. 71 ). 

Here we attempt a more intuitive derivation. For a very small source, we can consider a thin 
source annulus with angle A</> (see Fig. [3]). For a point lens, this thin annulus will be mapped into 
two annuli, one inside the Einstein ring and one outside. 

The area of the source annulus is given by the product of the radial width and the tangential 
length dr s x r s A<j>. Similarly, each image area is dr x rA<j>, and the magnification is given by 



dr x rAcf) r dr 



dr s x r s A(j) r s dr s 
For two images given in Equation ©, one finds 



(12) 



(r s + Vr s 2 + 4) 2 (r s - ^r a 2 + 4) 2 

4r s Vf s 2 + 4 4r s vV s + 4 

The magnification of the '+' image is positive, while the '— ' image is negative. The former image is 
said to have positive parity while the latter is negative (for the geometric meaning, see Fig. [3). The 
total magnification is given by 

/x=|/i + | + | M -|= r ' r ^- J , (14) 
r s Vr s z + 4 

and the difference is identically equal to unity 

|/i + |-|/x-| = l. (15) 

We make some remarks about the total magnification and image separations: 

1. When r s = 1, fj, = 3/\/5 w 1.342, corresponding to about 0.319 magnitude. Such a differ- 
ence is easily observable (For bright stars, the accuracy of photometry can reach a few milli- 
magnitudes.), and so the area occupied by the Einstein ring is usually taken as the lensing 
"cross-section." 



3 Imagine looking at a piece of white paper with a magnifying glass, the area is magnified, but the brightness of the paper 
is the same. 
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Fig. 3 Images of a thin annulus from r s to r s + dr s by a point lens on the plane of the sky. 
The dashed line is the Einstein ring. A<j) is the angle subtended by the thin annulus. Notice 
that the positive image (outside the Einstein ring) is clockwise while the negative image 
(inside the Einstein ring) is counter-clockwise. 



2. When r s — > oo, — s- ?' s 4 , fj, — > 1 + 2r s 4 . The angular image separation is given by 
A0 = (r s +2r s - 1 )0 E . 

3. High magnification events occur when r s 0. The asymptotic behaviors are [i — > ?' s _1 (l + 
3r s 2 /8), A6> -» (2 + r s 2 /4)9 E , and dr/dr s -> 1/2. The last relation implies that, at high 
magnification, the image is compressed by a factor of two in the radial direction (see Fig.|3}. 



2.5 Light curve and microlensing degeneracy 

Given a source trajectory, we can easily describe the standard light curve with a few parameters 
which suffers from the microlensing degeneracy. 

2.5.1 Standard light curve 

For convenience, we put the lens at the origin, and let the source move across the line of sight along 
the x-axis (see Fig. 0). The impact parameter in units of the Einstein radius is labeled as uq. For 
convenience, we define the Einstein radius crossing time (or 'timescale') as 

, r E 9 E r E _ v t 

*E = — = , o E = —, ^ rc i = — , (16) 

ft Mrol Dd D d 

where t> t is the transverse velocity (on the lens plane) and /j, rc \ is the relative lens-source proper 
motion. Substituting the expression for the Einstein radius into Equation (fTTT i. we find that 



t E « 19 day, LHTM (*) (M) ( * , ) 
y V D s \ D s J\8kpcJ \0.3MqJ V200kms- 1 
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Fig. 4 Illustration of the lens position and source trajectory. The dimensionless impact 
parameter is uq, (x s , y s ) is the dimensionless source position along the trajectory, and r s 
is the distance between the lens and source. 




If the closest approach is achieved at time t = to, then the dimensionless coordinates are x s 
(t — £o)Ae> Us = u o an d the magnification as a function of time is given by 

Two examples of light curves are shown in the right panel of Fig.Q]for Uq = 0.1 and 0.3. 

To model an observed light curve, three parameters are present in Equation ( fT8l : to, £e, and 
uq. In practice, we need two additional parameters, too, the baseline magnitude, and J^, a blending 
parameter. .F s characterizes the fraction of light contributed by the lensed source; in crowded stellar 
fields, each observed 'star' may be a co mposite of the lensed star, other unrelated s tars within the 
seeing disk and the lens if it is luminous (lAlcock et al. l200Tbl iKozIowski et al.l2007l) . Blending will 
lower the observed magnification and in general F a depends on the wavelength, and so each filter 
requires a separate F s parameter. 

Unfortunately, we can see from Equation ( TT8l that there is only one physical parameter (£e) in 
the model that relates to the lens' properties. Since depends on the lens mass, distances to the lens 
and source, and the transverse velocity /i rc i, from an observed light curve well fitted by the standard 
model, one cannot uniquely infer the lens distance and mass; this is the so-called microlensing 
degeneracy. However, given a lens mass function and some kinematic model of the Milky Way, we 
can infer the lens mass statistically. 



2.6 Non-standard light curves 

The standard model assumes the lensed source is point-like, both the lens and source are single 
and all the motions are linear. The majority (~ 90%) of microlensing events are well described by 
this simple model. However, about 10% of the light curves are non-standard (exotic), due to the 
breakdown of one (or more) of the assumptions. We briefly list these possibilities below. These non- 
standard microlensing events allow us to derive extra constraints, and partially lift the microlensing 
degeneracy. Because of this, they play a role far greater than their numbers suggest. 
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(1) binar y lens events: The lens may be in a binary or even a multiple system dMao & Paczriskil 
Il99lh . The light curves for a binary or multiple lensing system can be much more diverse 



(see Sec.|2.7l). They offer an exciting way to discover extrasolar plane ts ( Mao & Paczriskfl 1991 



iGould & Loebll 19921: iBennett & Rhie1ll996t iGriest & Safizadehlll998t iRattenburv etaU2002 , 



for more see Sec. 13. 51 ). 

(2) binary source events: The source is in a binary system. In this case, the light curve will be a 
simple, linear supe rposition of the two sources (when the orbital motion can be neglected, see 
lGriest&Hull992l) . 

(3) finite source size events: The finite size of the lensed star cannot be neglected. This occurs when 
the impact parameter uq is comparable to the stellar radius normalized to the Einstein radius, 
uo ~ r^/rs- In this case, the light cu rve is significantly modified by the finite source size 
effect dWitt & Mao 111994 lGould1l994l) . The finite source size effect is most important for high 
magnification events. 

(4) parallax/" xallarap" events: The standard light curve assumes all the motions are linear. 
However, the source and/or the lens may be in a binary system, furthermore, the Earth revolves 
around the Sun. All these motions induce a ccelerations. The effect due to Earth's motion aroun d 
the Sun is usually called "parallax" (e.g. IGould Ill992t ISmithl 120021: iPoindexter et al. Il2005l) . 
while t hat due to binary motion in the sourc e plane is called "xallarap" ("parallax" spelt back- 



wards, lBennettlfl998t lAlcock et al. I l2001al) . Parallax or "xallarap" events usually have long 



timescales. For a typical microlensing event with timescale £e ~ 20 day, the parallax effect due 
to the Earth rotation around the Sun is often undetectable (unless the photometric accuracy of 
the light curve is very high or the lens is very close). 
(5) Repeating events: Micro lensing can "repeat", in particular if the lens is a wide binary 
Stefano & Mao II 19961) or the source is a wide binary. In such cases, microlensing may man- 
ifest as two well-separated peaks, i.e., as a "repeating" event. A few percent of mi crolensing 
events are predicted to repeat, consistent with the observations dSkowron et al. 1 120091) . 

Notice that several violations may occur at the sa me time, which i n some cases allow the microlens- 
ing d egeneracy to be completely removed (e.g. lAn et alj 120021: iDong et all I2009b1 : iGaudi et al. I 
120081) . 

In particular, wh en both finite source size and parallax effects are seen, the lens mass can be 
determined uniquely dGould II 19921) . Microlensing parallax measures the projected Einstein radius in 
the observer plane (in units of the Earth-Sun separation, AU): Pb — J"e x D s /Da s , or equivalently, 
the parallax 7Te = 1 AU/rE- The finite source size effect, on the other hand, measures the ratio of the 
angular stellar radius 0* to the angular Einstein radius 8 y. Since the ang ular size of the lensed star 
can be measured from its color and apparent magnitude dYoo et alJl2004l) . we can derive the angular 
Einstein radius 6*e- It is straightforward to combine Eqs. dTTT ) and (O to obtain 

<19) 

Notice that the determination is independent of all the distances. 

Equation j\9[ is especially transparent in the natural formalism advocated by IGould I d2000l) . 
which provides a way to connect quantities measurable in microlensing (ttb, #e, ^e) with other phys- 
ical quantities: 

TTrcl = 7T E #E, Ihcl = — , (20) 
IB 

where tte = 1/^e, 7r r ci = V-Dd — l/D s is the lens-source relative parallax and yu rc i is the relative 
proper motion (also used in equation [T6b. For example, if the lens and distance can be measured 
(using other means), then 7r re i can be determined. In addition, if the lens and source motions can be 
measured, then we can determine the relative proper motions (/J. re i)- 
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2.7 iV-point lens gravitational microlensing 

It is straightforward to derive the (dimensionless) lens equation for TV-point lenses. We can first cast 
Equation ^ in vector form and then rearrange 

r s = r- tKoT. (21) 

| r |2 

The above expression implicitly assumes that the lens is at the origin, and all the lengths have been 
normalized to the Einstein radius corresponding to its mass (or equivalently, the lens mass has been 
assumed to be unity). 

Let us consider the general case where we have iV-point lenses, at r p. = (xk , Uk ) with mass 
k = 1, • • •, N. We normalize all the lengths with the Einstein radius corresponding to the total mass, 
M = J2k=i The generalized lens equation then reads 

N 

r s = r-> ?7i fc - r^, m k = —, (22) 

where 5Zfc=i m k = 1- If there is only one lens (mi = 1) and the lens is at the origin, then Equation 



j reverts to the single lens equation (12 It . 

Two-dimensional vectors and complex numbers are closely related, I Witt Id 19901) first demon- 
strated that the above equation can be cast into a complex form by direct substitutions of the vectors 
by complex numbers: 

N N 

Ez — Zk \ mk 
mk-. r— —r=z-]>- — , (23) 
fc=1 {z-Zk){z-Zk) t^\ z ~ Zk 

where z = x + y i, Zk — Xk + Uk i, and z s = x s + y s i (i is the imaginary unit). 

We can take the conjugate of Equation (|23l and obtain an expression for z. Substituting this back 
into Equation d23l ), we obtain a complex polynomial of degree TV 2 + 1. This immediately shows that 

1) even a binary lens equation cannot be solved analytically since it is a fifth-order polynomial, and 

2) the maximum number of images cannot exceed iV 2 + 1. In fact, the upper limit is 5(7V — 1) (see 
Sec. 13. 71 ). which indicates that most solutions of the complex polynomial are not real images of the 
lens equation. 

The magnification is related to the determinant of the Jacobian of the mapping from the source 
plane to the lens plane: (x s , y s ) —> (x, y). In the complex form, this is dWitt 111 9901) : 

1 j = d(x s ,y s ) = 1 dz s dz s 
J' d(x, y) dz dz ' 

Notice that the Jacobian can be equal to zero implying a (point) source will be infinitely magnified. 
The image positions satisfying J = form continuous critical curves, which are mapped into caus- 
tics in the source plane. Of course, stars are not point-like, they have finite sizes. The finite source 
size of a star smooths out the singularity. As a result, the magnification remains finite. 
For iV-point lenses, from the complex lens equation d23l , we have 



q N N 

Oz s x - TO fe _ _ | ^ m k 



k 

It follows that the critical curves are given by 



dz f^iz-Zk) 2 ' \^-' i {z - ZkY 



(25) 
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t/t E 

Fig. 5 Left panel: Caustics (red curves) and critical curves (green curves) for a binary 
lens. The lenses (indicated by two '+' signs) are located at z\ = (—0.665,0) and Z2 = 
(0.035, 0) with mass mi = 0.95 and TO2 = 1 — mi = 0.05 respectively. The black line 
shows the trajectory for three source sizes, r s /r^ = 0, 0.05, 0.3, indicated by the cyan and 
blue circles and a dot (for a point source). The trajectory starts at (—2,-1) with a slope of 
0.7. Right panel: Corresponding light curves for the three source sizes along the trajectory 
in the left panel. The time is normalized to the Einstein radius crossing time, t®, and t = 
corresponds to the starting position. Notice that as the source size increases, the lensing 
magnifications around the peaks usually decrease. 



The sum in the above equation must be on a unit circle, and the solution can be cast in a parametric 
form 

N 

EA^* (27) 

where < $ < 2tt is a parameter. The above equation is a complex polynomial of degree 2N with 
respect to z. For each $, there are at most 2N distinct solutions. As we vary <E> continuously, the 
solutions trace out at most 2 TV continuous critical curves (critical curves of different solutions can 
smoothly join with each other). In practice, we can solve the equation for one <E> value, and then 
use the Newton-Raphson method to find the solutions for other values of $. For a poin t source, the 
comp lex polyn omial can be easily solved numerically (e.g. using the zroots routine in lPress et al. I 
1992. see also ISkowron & Gould II20I2I However, for a source with finite size, the existence of 
singularities makes the integration time-consuming (see Sec.|4]i. The right panel in Fig. [5] shows 
the light curves for three source sizes along the trajectory indicated by the straight black line. As 
the source size increases, the lensing magnifications around the peaks usually decrease and become 
broader. Furthermore, the number of peaks may differ for different source sizes. 

2.8 Optical depth and event rates 

So far we have derived the lens equation and light curve for microlensing by a single star. In reality, 
hundreds of millions of stars are monitored, and rj 2000 unique microlensing events are discovered 
each year (OGLE-IV alone identifies about 1500 microlensing events each year). Clearly we need 
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some statistical quantities to describe microlensing experiments. For this, we need two key concepts: 
optical depth and event rate. 



2.8.1 Optical depth 

The optical depth (lensing probability) is the probability that a given source falls into the Einstein 
radius of any lensing star along the line of sight. Thus the optical depth can be expressed as 

n(D d ) (7rr E 2 ) dD A , (28) 

which is an integral of the product of the number density of lenses, the lensing cross-section (= 
7ITe 2 ) and the differential path (dD d ). 

Alternatively, the optical depth can be viewed as the fraction of sky covered by the angular areas 
of all the lenses, which yields another expression 

i r D ° 

t = — [n{D d )^D d 2 dD d ] (ir8 E 2 ) , (29) 
4?r Jo 

where the term in the [ ] gives the numbers of lenses in a spherical shell with radius D d to D d + dD d , 
TrOft 2 is the angular area covered by a single lens, and the term in the denominator is the total solid 
angle over all the sky (Air). 

If all the lenses have the same mass M, then n(D d ) = p(D d )/M, irr-E 2 oc M, and the lens 
mass drops out in n(D d )irr-E 2 . Therefore the optical depth depends on the total mass density along 
the line of sight, but not on the mass function. 

Let us consider a simple model where the density is constant along the line of sight, p{D d ) = po. 
Integrating along the line of sight one finds 

_ 2ttG Po a _ 1 Gpq4ttD s 3 /3 _ 1 GM(< D B ) _ V 2 
T 3c 2 s ~2c 2 D s _ 2c 2 D s ~2c 2 ' ( ' 

where M(< D s ) is the total mass enclosed within the sphere of radius D s and the circular velocity 
is given by V 2 = GM(< D s )/D s . 

For the Milky Way, V ~ 220 km s _1 , and r rj 2.6 x 10~ 7 . The low optical depth means millions 
of stars have to be monitored to have a realistic yield of microlensing events, and thus one needs to 
observe dense stellar f i elds, which i n turn means accurate crowded field photometry is essential (see 
lAlard & Luptodll998l:lAlardl2000tlWozni^l2000ll2008l and references therein). 



2.9 Event rate 

The optical depth indicates the probability of a given star that is within the Einstein radii of the lenses 
at any given instant. As such, the optical depth is a static concept. We are obviously interested in 
knowing the event rate (a dynamic concept), i.e., the number of (new) microlensing events per unit 
time for a given number of monitored stars, TV*. 

To calculate the event rate, it is easier to imagine the lenses are moving in a background of static 
stellar sources. For simplicity, let us assume all the lenses move with the same velocity of v t . The 
new area swept out by each lens in the time interval dt is equal to the product of the diameter of the 
Einstein ring and the distance traveled v t dt, dA = 2?'e x v t dt = 2r-E 2 dt/tE- The probability of a 
source becoming a new microlensing event is given by 

dr = n(D d )dAdD d = n(D d ) ) dtdD d (31) 
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The total number of new events is N+dr, and thus the event rate is given by 

- w .r BW f-Lww«w?. (32) 



dt Jq / T J t~E 

If, for simplicity, we assume that all the Einstein radius crossing times are identical, then we have 



7T ts 

We make several remarks about the event rate: 



(33) 



(1) If we take = 19 day (roughly equal to the median of the observed timescales), then we have 
r « ™tL = 1200yr-^^ ( " . (34) 



7T t E 10 8 lO- 6 V 19 day 

For OGLE-III, about 2 x 10 s stars are monitored, so the total number of lensing events we 
expect per year is T ~ 2400 if r ~ 10~ 6 , which is within a factor of four of the observed rate 
(indicating the detection efficiency may be of the order of 30%). 

(2) While the optical depth does not depend on the mass function, the event rate does because of 
tE(oc M 1 ' 2 ) in the denominator of Equation d33b . The timescale distribution can thus be used 
to probe the kinematics and mass function of lenses in the Milky Way. 

(3) The lenses and sources have velocity distributions; one must account for them when real- 
istic event rates are needed. Furthermore, the source distance is unknown, an d so in gen- 
eral we need to average o ver the source distance (for example calculations, see lGries"3 1 1 99 lfc 
iKirasa & Paczvriski II 19941) . 



3 APPLICATIONS OF GRAVITATIONAL MICROLENSING 

As we can see, the theory of gravitational microlensing is relatively simple. When IPacz vnskil d 1986b 
first proposed microlensing, the astronomical community was very skeptical whether microlensing 
events could be differentiated from other types of variable stars (for an example of skepticism, see 
lMadl2008bh . As proven again and again in the field, such pessimism is often unwarranted: with the 
rapid increase in the discovery rate of microlensing events (from a few events in the early years 
to about 2000 events annually), even small probability events can be discovered. The b e st exa mple 
may be the s e eming ly crazy idea of terrestrial parallax propos ed by Hardy & Walkerl d 1995b and 



iHoIz & Waldl dl996l) . which has subsequently been observed by iGould et al. ( 2009 ). Since the dis 



covery of the first microlensing events in 1993, more than 10,000 events (mostly in real-time) have 
been discovered. In the process, tens of TBytes of data have been accumulated. This tremendous 
database is a treasure trove for exploring diverse astrophysical applications. 



3.1 MACHOs in the Galactic halo? 

Microlensing was originally proposed to detect MACHOs in the Galactic halo jPaczvriskill 19861) . 
Earlier results by the MA CHO collaboration suggest that a substantial fraction (« 20%) of the 
halo may be in MACHOs dAlcock e t al. 2000) based on 13-17 events toward the Large Magellanic 
Cloud (LMC) in 5.7 years of data. The number of events qu oted here varies somewhat depending 
on the selection criteria , which is a debated point (see, e.g.. iBelokurov et al] 120041 ; iBennettl 12005 : 
iGriest & Thomas 112005b [Evans & Belokurovll2007l) . 

In any case, the result tu rns out to be controversial, as these numbers are not confirmed by the 
EROS dTisserand et al.ll2007l) or OGLE collaborations dWvrzvkowski et al.ll20lTl) . A recent analy- 
sis of the OGLE data toward the LMC concluded only <; 2% of the halo could be in MACHOs 
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(tWyrzykowski et al1l201 ll) - all events can be explained by lensing by normal stars. In fact, at least 
two of the MACHO events tu rn out to be by stars in the thick disk of the Milky Way dDrake et alj 
l2004lKallivavalil et alj2003) . The lack of MACHOs is entirely consistent with the baryonic content 
deter mined from many oth er astrophysical observations (e.g. from microwave background observa- 
tions, iKomatsuetaDIIoill). 



3.2 Galactic Structure 

Microlensing data offer a number of ways to study the Galactic structure, which is still somewhat 
under-explored, and more work is needed. 



3.2.1 Color-Magnitude Diagrams 

Microlensing surveys yield many high-quality color-magnitude diagrams of stellar populations over 
hundreds of square degrees. Such data contain much i nformation abou t stellar populations in the 
bulge. Despite of some promising earlier attempts (e.g. lNg et al. Ill996h . not much work has been 
done since due to the difficulty of disentangling the dust extinction, spatial distributions and com- 
plexities in the stellar populations (see the next subsections). It is nevertheless worthwhile to perform 
further analyses in this area. 



3.2.2 Red Clump Giants As Standard Candles 

Red clump giants are metal-rich, core-helium bur ning horizontal branch sta rs. They have ap 



proximately constant luminosity (standard candles) dPaczvhski & Stanek 



1998b and relatively lit 



tie dependence on the metallicities, especially in the /-band (see IUdalskilll998at IZhao et aT 
iNataf & Udalskfeoi ll and references therein). As such, they can b e used to determine the distances 
to the Galactic center and the Magellanic Clouds dUdalskd 1 1 9 9 8bl) . This complements the determi - 
nation of other standard candles such as RR Lyrae stars and Cepheids dGroenewegen et aT1l2008l) . 
The distance to the Galactic center determined with these methods is typically very close to 8 kpc, 
with a combined systematic and statistical error of about 5%. 

The observed width of red clump giants in a given line-of-sight toward the Galactic center is 
larger than their intrinsic scatters (w 0.2 mag) because of the finite depth of the Galactic bar. A 
careful analysis of the counts of red clump giants can be used to determine the geometric parameters 
of the bar. For example, for 44 fields from the OGLE-II data, the three a xial lengths are found to be 
close to 10:3.5: 2.6 with a bar angle of 24-27° dRattenburv et al. Il2007al) . largely consistent with the 
earlier study bv lStanek et al.l dl997l) . Extra care, h owever, needs to b e taken for fields off the plane 



due to the presence of two red clump pop ulations (INataf et al. 2010l). which suggests an X-shaped 



structure in the Milky Way center (see also iMc William & ZoccarfeOlOHSaito et alj|201 lb . 

We note that the OGLE-III and OGLE-IV campaigns cover much larger areas on the sky, which 
can be used to constrain not only the bar but perhaps also spiral arm structures. Investigations taking 
this approach are already under way (Nataf et al. 2012, in preparation). 



3.2.3 Extinction Maps 

Red clump giants can also be used to infer extinctions tow ard the Galacti c Center. The first applica- 
tion of this m e thod t o the OGLE-I data was performed by IStanekl d 19961) and later to the OGLE-II 
data bv lSumil d2004 . The extinction ma ps are approximately consistent with those obtained with 
other methods (e.g. lSchlegel et ai]|l998l) . One interesting conclusion from these stu dies is that the 
extinction law is almost always anomalous toward the Galactic center dUdalskill2003l) . 

A very rece nt application of the m ethod to the VISTA Variables in the Via Lactea (VVV) survey 
can be found in lGonzalez et alj d2012l) . 
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3.2.4 Proper Motions 

Nearly two decades of time series of observations for many different bulge fields offer a way to 
determine stellar proper motions. This has been done for th e O GLE-II catalog for about f ive million 
stars with an accuracy of ~ mas yr _1 (see lSumi et al.l2004l and lR~attenburv & Mao l2008l for a study 
of selected high-proper motion stars). 

One particularly interesting issue concerns the anisotropy in the bulge kinematics. The proper 
motion dispersions in the longitudinal and latitudinal directions have ratios close to w 0.9; these 
values are similar to those independently measured with the Hubble Space Telescope in a number 
of fields but for a smaller number of stars dKozlowski et al.ll2006t Idarkson et al.ll2008l) . A compar- 
ison with theoretical pred ictions seems to indicate that the theoretical models are too anisotropic 
dRattenburv et al. 2007bl).T his effect is also seen in more recent modelling with the Schwarzschild 



method (IWang et al. 112012b . The reason for this discrepancy is not completely clear. One issue worth 



mentioning is that bulge observations are often contaminated by disk stars, a point emphasized by 
IVieira et alJd2007l) . 



3.2.5 Optical Depths and Galactic bar 



In Sec. 12. 81 we showed that the optical depth provides an independent probe of the density distri- 
bution in the bulge. Earlier determinations of the optical dep th in the Baade Window have provided 
indep endent evidence that the Galactic bulge contains a bar dKiraga & Paczvriski l ll994tlZhaoeTai1 
119961) . While we have discovered about ten thousand microlensing events over the last two decades, 
only a small fraction has been used for statistical analyses of optical depths. 

More recent deter minations of microlensing optical depths have been performed by 
lHamadache et al.l (12006b using 120 red cl ump giants from the E ROS collaboration, 60 red clump 
giants f rom the MACHO collaboration dPopowski et al.l 120051) . 28 events for the MOA collab- 
oration ( Surni et all 120031) and 32 high signal-to-noise ratio events for the OGLE collaboration 



(ISumi et al. 2006b. The measured optical depths are largely consistent with theoretical predictions 
(e.g. JWood & Maoll2005l: iRvu et alJ2008t iKerins et alJl2009b . 



The reason that the measurements have only been performed for red clump giants is that they 
are bright and are supp osed to suffer less from the blending in the crowded bulge field (in fact, this is 
not rigorously true, see lSumi et al 120061) . If we can overcome this and the issue of human resources 
(see the next subsection), combined with the number counts of red clump giants, we can in principle 
constrain the bulge density distribution much better using non-parametric models. For example, we 
can use non-parametric models rather than the simple parametric models often used in the l iterature 
(still based on very poor-resolution [sa 7°] maps from DIRBE on COBE. lDwek et afll 1995b . 



3.2.6 Timescale Distributions and Mass Functions 

As we have seen in Sec. 12. 91 the timescale d istribution of event ra tes carries information about the 
mass function of lenses. An early study by lHan & Gouldl (1 1995b of 51 events concluded that for 
a power-law mass function, n(M)dM ~ M~@ dM, the preferred slope (3 w 2.1, with a lower 
m ass cutoff o f abou t 0.04M Q . Their slope is close to the Salpeter value {(3 = 2.35). The study 
bv lZhao et al. (1996) also concluded that the data are inconsistent with a large population of brown 
dwarfs in the bulge. A more recent analysis bv lCalchi Novati et al. I d2008l) modeled the mass function 
andfound/3 w 1. 7±0. 5. The large error bar is again due to the small number of events used (~ 100). 

The current data base contains roughly two orders of magnitude more events - it would be very 
interesting to explore what we can learn with the full data set we now possess. However, there 
are at least two difficulties that will need to overcome. One is human resources: extensive Monte 
Carlo simulations must be performed to determine the completeness of surveys - this will be a time- 
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consuming exercise. Secondly, the effect of blending of background stars needs to be taken into 
account (s ee Sec . 12. 5. lb . This is not so straightforward since the fields have different degrees of 
crowding dSmith et alj|2007l) . Nevertheless, it will be valuable to perform an analysis of the optical 
depth maps and timescale distributions using all the data set to understand the density map and 
mass functions of the bulge and mass of the disk. Notice that the microlensing sample is largely a 
mass-selected sample, rather than a light-selected sample as in most other studies. 

This exercise may be particularly interesting in light of evidence for systematic variatio ns in 
the initial mass functions i n elliptical galaxi es from both dynamics dCappellari et al.ll2012l) and 
strong gravitational lensing dTreu et al.l 2010h . A careful study of the mass functions in the bulge 
and disk from microlensing offers an important independent check on these conclusions by investi- 
gating whether the Galactic bulge follows similar trends. 



3.3 Stellar atmosphere and bulge formation 

High-magnification events offer great targets-of-opp ortunity to obtain sp ectra with high signal-to- 
noise ratio to study stellar atmospheres in bulge stars dLennon et al.l 19960 . More recent observations 
using 8m class telescopes allow determinations of a number of stellar parameters such a s surface 
gravity, metallicity and ages (e. g.lLennon etdll996tlThurl et al.ll2006UCohen et al.1l2008l) . 

The most recent studies bv lBensbv et al.l d2010il201l1) of 38 stars indicate that for dwarf stars, 
there may be two populations of stars. The metal-poor population seems to be predominantly old 
(age ~ 12 Gyr) while the metal-rich population appears to have a bi-modal distribution: one is old 
(~ 12 Gyr) while the other one is intermediate age (3-4 Gyr). The existence of int ermediate age 
stars seems to be in conflict with broad-band color and spectroscopy of giants (e.g. IZoccali et alj 
120081) . One way to resolve this confl ict is that the age estimate may be incorrect due to enhanced 
helium enrichment dNataf et al.ll20Tlb . While more microlensed stars are desirable, it demonstrates 
the power of microlensing as a natural telescope, similar to clusters of galaxies as a natural telescope 
to study very high-redshift galaxies. 

The surface brightness of a star is not uniform, instead its limb usually appears darker ("limb- 
darkening"). Limb-darkening profiles have been measured for the Sun and a few other stars. The 
sharp magnification gradient in high-magnification or caustic-crossing events allows us to study the 
limb-darkening profile as the source moves across the line of sight or caustics. Notice that limb- 
darkening profiles can not only be studied in br oad-band but also in spectral lines such as Ha if one 
has time-resolved spectra during microlensing dThurl et al.ll2.Q06h . 

Approximately 10 G and K giants had their limb-dar kening profiles measured with microlensing 
dCassan et al.1 l2006t iThurlet al.l 120061; IZub et al.l 1201 lb . These s tudies suggest that "the classical 
laws are too restrictive to fit well the microlensing observations" dCassan et alj|2006l) . and radiative 
transfer models will need to be improved. 



3.4 Mass Determinations of Isolated Stars: from Brown Dwarfs to Stellar Mass Black Holes 

As shown in Sec. 12.5.11 for standard microlensing events, the lens mass cannot be determined 
uniquely. However, for exotic events, partial or complete removal of this degeneracy is possible (see 
equation [T9l. Microlensing thus provides an important new way to determine the mass of isolated 
stars. 

All stellar black hole candidat es in the Milky Way are in bina ry systems and have been discov- 
ered through X-ray emissions (see iRemillard & McClintockll2006l for a review). Their masses range 
from 5Mq to 30M©. Microlensing provides an independent way to study isolated black holes. The 
principle is very simple: Typical microlenses have masses of about 0.3A/©, black holes are ^ 10 
more massive, and so events due to stellar mass black holes should be a factor of a few longer. These 
events thus have a much larger chance of exhibiting parallax signatures, which can be used to de- 
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termine the lens as a function of distance. Combined with a mass density and kinematic models of 
the Milky Way, we can constrain the lens mass. If it is higher than a few solar masses and it is dark 
(as can be inferred from the light curve from the blending parameter ,F S , see Sec. 12.5.1k then it is a 
potential stellar mass black hole candidate. 

Hal f a dozen or so such candidates have been i dentified by dMao et al .120021: iBennett et al.l2002l 
see also lAgol et al.l 120021: iPoindexter et aL~ll2Q05h . A searc h for the X-ray emission in one of th e 
candidates, MACHO-96-BLG-5, yielded only an upper limit dMaeda et al.l2005l:lNucita et al.l2006l) . 
consistent with the system being a truly isolated black hole. 

An ambitious 3-year HST survey (192 orbits) is under way by a team led by K. Sahu to detect 
microlensing events caused by non-luminous isolated black holes and other stellar remnants. It will 
be very interesting to see the results from this survey (ISahu et al.ll2.Q12b . 

Other direct mass measurements have also been performed, especially for h ighly magnified mi- 
crolensing events. The determined ma sses range from brown dwarf c andidates ( Smith et alj|2003bt : 
iGould et"aill2009l iHwang et alj|2010l) to more normal stars (see, e.g. jBatista et al. 20091 and refer- 
ences therein). 



3.5 Extrasolar Planets 



Undoubtedly the highlight of gravitational microlensing in the last decade has been the discov- 
ery of extrasolar planets. At the time of writing, about 20 microlensing extraso l ar planets have 



been 


discovered, including 16 published dBondetal. 20041 Udalsk et al.l l2005t Beaulieu et al. 


2006 


; IGould et al. 1 20061; Bennett et al. 1 2008t Gaudi et al. 2008F Dong et all 2009b; Janczak et al. 


2010 


; Isumi et al. 20101 iMivake et al. 201 1 ; Batista et all 201 1 ; Muraki et al. 201 ll Yee et al.l 2012; 


Bachelet et al.l 2012; iBennett et al. I2OI2I). Of these, 11 are high-magnification events, demonstrat- 



ing they are excellent c andidates for hunting planets, as pointed out by dGriest & Safizadeh"lll998t ; 
IRattenburv et al. Il2002l) . While this is a small fraction of the 800 extrasolar planets discovered so 
faiQ, they occupy a distinct part of the parameter space which would be difficult to access with other 
methods (see below). 



The method itself wa s proposed more t h an 20 years ago dMao & Paczhskilll99ll;lGould & Loebl 
119921) . In the abstract of iMao & PacznskH dl99ll) . it was optimistically claimed that 'A massive 
search for microlensing of the Galactic bulge stars may lead to a discovery of the first extrasolar plan- 
etary systems." Paczyriski, however, mentioned the idea as "science fiction" at the 1991 Hamburg 
gravitational lensing conference. In reality, it took mo re than a decade of microlensing observations 
for the first convincing extrasolar planet to be found dBond et aT"ll2004l). with heroic e ffort s in be- 
tween . Much of the theory and observations have been reviewed by IRattenburv! d2006l) and iGaudTI 
d2010h : we refer the readers to those papers for further details. 

Fig. [6] shows an example of an extrasolar planet discovered by microlensing. The extrasolar 
planet has a mass of around 5.5M0, manifested as a secondary bump on the declining wing of the 
light curve, lasting for about one day. This illustrates that to find an extrasolar planet, the dense 
sampling of light curves plays a critical role. 

The first two-planet system discovered bv lGaudi et al. I d2008l) through monitoring of high mag- 
nification events is shown in Fig. [7] The light curve in this case is much more d ramatic due to the 
complex caustics involved (the top left inset) and orbital motion in the system (see lPennv et al.ll2~01 ll 
for detailed predictions). The mass and separations of the two planets are very much like those for 
the Saturn and Jupiter in our solar system. Such multiple planet systems constitute about 10% of 
the extrasolar planets discovered through microlensin g. Many multiple planetary systems have b een 
discovered in radial velocity and transit surveys (e.g.. IWright et alj|2009l iFabrvckv et alj|2012l and 
references therein). In radial velocity surveys, at least 28% of known planetary systems appear to 
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Fig. 6 A super-Earth (« 5.5M(q) discovered by microlensing. Different symbols indicate 
data from different observatories; the solid line shows the best-fit model(s). The two insets 
show the OGLE data alone and the planetary bump lasting for about a day. Reprinted by 
permission from Macmillan Publishers Ltd: Nature 439: 437-440, ©2006. 



contain multiple planets (IWrightetal.ll2009h . However, these two fractions cannot be easily com- 
pared since the planets discovered are at very different separations from the host stars. Furthermore 
microlensing is only sensitive to planets within a narrow range of the Einstein radius. A more de- 
tailed comparison is needed to address this issue. 

Table [T]lists all the published extrasolar planets discovered b y microlensing at the time of writ- 
ing. This is an expanded version of Table 3 in lMivake et al.ld201 lb . The microlensing planets occupy 
a distinct part in the plane of mass vs. separation normalized by the snow line. This is illustrated in 
Fig- El It is seen that microlensing planets reside mostly outside the snow line where their equilib- 
rium temperatures are low, about 100 K (see also, e.g., Fig. 3 in lBeaulieu et al. Il2008l) . Beyond the 
snow line, the density of solid particles in the disk increases rather abruptly, the planetesimal cores 
can form faster and the formation of gas giants becomes easier. The equilibrium temperatures are 
very different from the hot planets found, for example, in radial velocity searches. 

Analyses of microlensing extrasolar planets has already yielded very interesting statistical re- 
sults on pl anet frequency. Su rprisingly, six high-magnification events appear to form a well-defined 
"sample" dGould et alj|2010h even though the observations by the followup teams were triggered 
by somewhat chaotic human involvement. The authors presented the first measurement of the 
planet frequency beyond the snow line, for the planet-to-star mass-ratio interval —4.5 < logg < 
—2, corresponding to the range of ice giants to gas giants. The frequency was found to follow 
d 2 N p \/d log q/d logs = (0.36 ± 0.15)de x~ 2 , where q is the m ass ratio and s is the separation. 
This is consistent with the extrapolation of ICumming et alj d2008[) to large separations. Their study 
also implies a first estimate of 1/6 for the frequency of solar-like systems. 

More recentlv. lCassan et al.ld2012l) reported a statistical analysis of microlensing data (gathered 
in 2002-07) and concluded that about 20% stars host Jupiter-mass planets (between 0.3 and 10 
Jupiter masses) while cool Neptunes (10 — 301/^) and super-Earths (5 — lOA/^) are even more 
common: their abundances per star are close to 60%. 
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Fig. 7 The first two-planet system discovered by iGaudi et al. I (120081) . Different symbols 
are data taken from different telescopes (12 in total). The numbers indicate different fea- 
tures produced by caustic crossings and orbital motions. From Science, Volume 319, no. 
5865, pp. 927-930. Reprinted with permission from AAAS. 



Somewhat controversially, the study by ISumi et al 1 J201 1) found a population of free-floating 
Jupiters, with 1.8 Jupiters per star on average. These planets are not bound to any parent stars, 
observationally they manifest as very short-time scale events (in contrast to very long events for 
stellar mass black hole candidates). Whether such a high frequency of free-floating Jupiter-mass 
planets can be produced in core accretion theory or gravitational instability theory is unclear. 

3.6 Numerical Methods in Modelling Binary and Multiple Light Curves 

The modeling of light curves can be naturally divided into two parts: The effi cient calculatio n of 
light curves and finding the best-fit models in a \ 2 sense. We refer the readers to lBenneti (1201 Ob for 
an excellent discussion of this topic. 

3.6.1 Calculation of light curves 

For a single lens, the image magnifications and positions are analytical. The light curve calculation 
is straightforward, even accounting for finite source size (IWitt&Maolll994l) . 

For a finite size source lensed by a binary or multiple system, its magnification can be calculated 
in three different regimes with increasing complexities. 
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Table 1 Parameters of exoplanets discovered by microlensing. 



Name 


Host Star Mass 


Distance 


Planet Mass 


Separation 
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a MOA-2()07-BLG-400Lb has two solutions due to a strong close/wide model degeneracy (see Sec . 13. 7) . 
b Details of the MOA-2008-BLG-310Lb parameters are discussed by Ijanczak eta?] feOlfjft and lSumi et a?] 4201Ct ). The 
en'or bars take into account the degeneracy. 
c 90% confidence limit. 

d The system has two solutions due to a close/wide separation degeneracy (see Sec . 13. 71 , thus two separations are given 
from a Bayesian analysis assuming the lens is a main sequence star. The close separation is slightly favored. 
Notes: LB: lens brightness; DL: detection of the lens. 

Referen ces: 1 .iBond et all <2004l): 2.lBennett et all hOOcl) : 3.IUdalsk et alj 42005b. 4.lDong et alj 12009 d) : 5. lGould et al" 
6. iBeaulieu et al. I fc006l) ~7. iGaudi et al. I 42008l); ~8. ISumi et alj feOlOl); 9 . iBennett et aTl <20O8ll: 10. iDong et ~ 



2009U) : l l lJanczaket^ koiOl) : 12 iMivake et alj 1201 
20121); 16. lBachelet et alj)2012l) ; 17. lBennett et"aiT)2012l) . 



When the source is far away from the caustics or any sharp magnification gradient, then the 
point source approximation can be used. In this case, the magnification can be readily found 
by solving the fifth-order or t he tenth-order poly nomials for binary or triple lens systems (e.g. 
using the z roots routine in Press et a for a more efficient implementation available 

publically see lSkowron & Gould 12012b . In practice, as the source moves along its trajectory, we 
can feed the image positions from the previous step as the initial guesses and use the Newton- 
Raphson (Simpson 1740) method to achieve rapid convergence to new solutions for the current 
source position. These can be used to deflate the polynomial to a lower-order one, which can 
then be solved much more easily. 

In the case of binary lensing, the deflated polynomial is usually either linear or quadratic and 
its solutions can then be found analytically. In this approach, we, at all times, keep all the five 
solutions to the polynomial (not all are true solutions of the le ns equation). This effi c ient m ethod 
was used in the first attempt to pre dict binary light curves Ijaroszvhski & Maol (1200 lb : it is 
similar to ISkowron & Gould I d2012l) in spirit. We note that for five-image configurations, the 
total magnification of positive-parity images must exceed that of negative-parity images by unity 
dWitt&Madll995h . This can be used as a test of numerical accuracy. 

When the source is moderately (one diameter or so) far away from the caustics, th e magnificatio n 
can be efficiently calculated using the hexadecapole approximation proposed bv lGouldl d2008l) . 
It has the advantage of being very efficient and at the sam e time providing an estimate of the 
approximate accuracy (see also lPeicha & Hevrovsky|2009l) . 

When the source is approximately within one diameter of the caustics, the magnification for a fi- 
nite source size calculation becomes complex and time-consuming due to the presence of singu- 
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semi-major axis/snow-line 



Fig. 8 Extrasolar planets in the plane of mass vs. separation (in units of the snow line, 
indicated by the vertical dashed green line). The snow line is taken to be at « 2.7Ml/M q 
AU, where Ml is the lens mass (see TableQ]). The red filled circles with error bars indicate 
planets found by microlensing. The open symbols show the degenerate close-separation 
solutions for MOA-2007-BLG-400Lb and MOA-201 l-BLG-293Lb (see Table 1). The 
black triangles and blue squares indicate the planets discovered by radial velocities and 
transits, respectively. The magenta and green triangles indicate the planets detected via 
direct imaging and timing, respectively. The non-microlensing exoplanet data were taken 
from The Extrasolar Planets Encyclopaedia (http://exoplanet.eu/). The planets in our solar 
system are indicated with initial letters (in cyan). 



larities. Many methods have been proposed, including (modified) rav shooting dRattenburv et a"L~1 
20021: iDong et al.ll2006l) . level contou ring (based on Stokes theorem. Ipould & Gaucherelll 1997c 



Dominik 2007b and grid integration ( Bennett & Rhie |[l996l : iBennett 20101) . For detailed com- 
parisons between these methods, see IBennett d2010h . It is clear that there is still room for im- 
provement. Furthermore, several independent modeling codes are desirable for cross-checks. 
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3.6.2 Finding the global minimum 

Once a binary or multiple light curve can be efficiently calculated, it still remains a highly non-trivial 
job to find the best-fit parameters, especially in the presence of degeneracy (see Sec. |3.7t . Starting 
from an initial guess of the parameters, many routines can be used to converge to a local minimum 
using, e.g., MINUIT0 or GSL routines^ . More recently the Monte Carlo Markov Chain method has 
gained popularity (e.g. lDong et al.ll2006t lBennettll2010l) . 

However, to find the global minimum, often a grid of initial guesses are used. This may be, how- 
ever, impractical once the dimension of the parameter space increases, especially for multiple lenses 
(for a nice discussion of the parameters in binary lensing including parallax and orbital motion, see 
ISkowron et al. l l20TH) . This may be particularly severe when th e next-generat ion experiments come 
online; this area deserves much more research in the future (see lBenne tt 2010 for the current state of 
affairs). 



3.7 Mathematics of Gravitational Microlensing 

Gravitational l ensing is mathemat ically rich and is related to singularity (catastrophe) theory (see the 
monograph bv lPetters et al.ll200~Th . 

In classical Keplerian potential theory, the two-body problem is analytically tractable. In gravi- 
tational lensing, as we have shown in Sec. 12. 71 it is not even possible to solve the binary lens equation 
analytically. Curiously, for binary lenses, there is still an analytical theorem that five -image configu- 
rations must have a magnification no smaller than 3 dWitt & Maolll995URhie 1997b . This was used 



to infer the presence of blending in gravitational micr olensing ( Witt & Maoll 1995b 



In binary lensing, there is a degeneracy found by iDominikl i 1 999b betwe en close and wide sep- 
aration binaries. This was later explored in much greater detai l by lAnl d2005l) . Planetary and binary 
lens light curves can also mimic each other dChoi et al.ll2012b . It is unclear whether some of these 
degeneracies can be generalized to multiple (N > 3) lenses. For parallax events, there are also 
degeneracies (ISmith et al. 112003 at iGould Ill995b . 

The number of critical curves for iV-point lenses cannot e xceed 2N (see Sec. | 2.7|>. The upper 



bound on the number of images for iV-point lenses is 5(N- 1) dRhiell2003l Khavinson & Neumann 



2001 



2006). It is unclear whether the se two linear scalings are related in a topological sense (see [Rhie 



120031: lAn& Evans 112006b 



4 FUTURE OF GRAVITATIONAL MICROLENSING 

The future of gravitational microlensing is bright. In terms of current experiments, OGLE-IV already 
has a field of view of 1.4 square degrees and is readily identifying about 1500 events in real-time per 
year. The rate can in principle be increased by ~ 40% by analyzing their archives. 

The 1.8m MOA-II telescope has a field of view of 2.2 square degrees and the MOA collabora- 
tion issues about 600 microlensing alerts each year. The MOA collaboration uses difference image 
analysis photometry to issue alerts, which can identify faint source stars that are undetectable when 
unmagnified. This is in contrast to the OGLE, which to date has only been triggered by stars iden- 
tified using their template image. The MOA strategy appears to increase the alert rate by roughly a 
factor of ~ 1.4 compared to OGLE. 

There has also been an influx of new telescopes involved in microlensing, including the WISE 
Observatory which is conducting an independent survey. In the future, more telescopes may engage 
in microlensing campaigns, including the SKYMAPPER. Also, Chinese astronomers are discussing 
possible monitoring from Antarctica (Dome-A) and Argentina where 2m class telescopes have been 



5 http://wwwasdoc.web.cern.ch/wwwasdoc/minuit/minmain.html 

6 http://www.gnu.org/software/gsl/ 
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proposed by the Chinese Academy of Sciences. The superior seeing and continuous time coverage 
for a few months per year at Dom e-A may be an adv antage for microlensing (although the seeing is 
degraded due to the high air mass, IWang et alJ2009l) . In fact, as the first step, the first of three 50 cm 

telescopes has already been installed at Dome-A in early 2012 and will perfor m pilot surveys . 

In terms of hunting for extrasolar planets, sever al White Papers dGould et al. I 120071 
iBennettet al~l 120071 15 omim kl l2008t iBeaulieu et al. ||2008b set out strategies with ambitious mile- 
stones for the next fifteen years. 

- Current extrasolar discovery with microlensing is a mixture of surveys and intense followup 
observations. However, the boundary between surveys and followup teams is already starting to 
blur, with survey teams for some fraction o f their time, engag ing in very dense monitorings of 
some fields, for example, by both the MOA (ISumi et al.ll201 1 ) and OGLE teams. 

- In the next five years, a network of three 1.6m telescopes will be built by Korean astronomers 
(Korean Microlensing Network, KMTNet). These telescopes will be sited in Chile, South Africa 
and Australia. Each is equipped with a 4 square degree field camera, and will be able to monitor 
16 square degrees and densely sample light curves (every 10 minutes or so). This will in principle 
get rid of the cumbersome division between survey and followup. As a result, the selection 
functions will be much better specified and thus statistical studies will be easier to perform. This 
is important as the number of microlensing extrasolar planets will increase significantly. 

- A microlensing telescope in space in the next 10-15 years has been proposed both in the 
US and Europe. WFIR ST is the top recommended space mission by the US Decadal Survey 
dBlandford et al.ll2010l) . The recently funded Euclid mission by ES A may also have a microlens- 
ing component (for detailed simulations, see lPennv et al. ||2012|) . Observing from space has sub- 
stantial advantages: smaller and more stable PSFs and continuous time coverage will allow us 
to search for lower-mass (Earth-mass) planets. Space missions will also allow us to uniquely 
determine the masses of many extrasolar planets. Combined with the stellar transit mission 
Kepler, space microlensing experiment(s) will provide a complete census of Earth-mass (and 
lower-mass) planets at virtually all separations, including free-floating ones. 

The theory of microlensing is well understood, although computationally there are still some 
challenging issues (see Sec. 13. 6b . For example, it is still time-consuming to calculate the light curves 
for finite-size sources since we need to integrate over the singularities of caustics. This is particularly 
important for the discovery of extrasolar planets when a source transits t he small caustics in duced 
by the planet(s). The problem becomes even worse for multiple planets (IGaudi et al. 112008b . How 
do we efficiently search the high dimensional parameter space? Are there hidden multiple planetary 
light curves in the database that are not yet identified due to their complex shapes? 

With a very healthy interplay between theory and observations, upgraded/new surveys in the 
near term and space satellites on the horizon, microlensing can expect another exciting decade in the 
future. 

Acknowledgements I thank Drs. Liang Cao, Lijun Gou, Andy Gould and Richard Long for many 
helpful discussions and criticisms on the review. 
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